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Abstract. We present new observations of the cluster- 
lens Abell 370: a deep HST/WFPC2 F675W image and 
ESO 3.6m spectroscopy of faint galaxies. These observa- 
tions shed new light on the statistical properties of faint 
lensed galaxies. In particular, we spectroscopically con- 
firm the multiple image nature of the B2-B3 gravita- 
tional pair (Kneib et al. 1993), and determine a redshift 
of z — 0.806 ± 0.002 which is in very good agreement 
with earlier predictions. A refined mass model of the clus- 
ter core (that includes cluster galaxy halos) is presented, 
based on a number of newly identified multiple images. 
Following Bezecourt et al. (1998a), we combine the new 
cluster mass model with a spectrophotometric prescrip- 
tion for galaxy evolution to predict the arclets number 
counts and redshift distribution in the HST image. In par- 
ticular, the ellipticity distribution of background sources 
is taken into account, in order to properly estimate the 
statistical number and redshift distribution of arclets. We 
show that the redshift distribution of arclets, and partic- 
ularly its high redshift tail can be used as a strong con- 
straint to disentangle different galaxy evolution scenarios. 
A hierarchical model which includes number density evo- 
lution is favored by our analysis. Finally, we compute the 
depletion curves for the faint galaxy number counts and 
discuss their wavelength dependence. 
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1. Introduction 

Thanks to gravitational amplification, galaxy clusters are 
a powerful tool to probe distant galaxies. Indeed, one of 
the most distant galaxy detected is a gravitational arc 
at z=4.92 in the cluster C11358+62 (Franx et al. 1997). 
Other very distant lensed sources (z >3) have been iden- 
tified behind C10939+4713 (Trager et al. 1997) and A2390 
(Pello et al. 1998a, Frye & Broadhurst 1998) leading to 
important results on the formation history and evolution 
of galaxies at large redshift. Determining the redshift of 
arcs and arclets is of great importance as it fixes the an- 
gular scales of the optical configuration, hence giving an 
absolute cluster mass estimate within the arc radius (e.g. 
Soucail et al. 1988, Mellier et al. 1991). But despite the 
cluster magnification, measuring arclets redshifts is a dif- 
ficult observational task due to their low surface bright- 
ness (Bezecourt & Soucail 1997, Ebbels et al. 1998), and 
the lack of strong spectral features in the optical domain 
for galaxies with 1 < z < 2.5. 

Accurate modeling of cluster potentials based on 
the analysis of multiple images and weak shear dis- 
tortions has shown that cluster mass distributions are 
best described by the sum of a smooth and large scale 
component (the cluster) and the contribution of clus- 
ter galaxy halos (Kneib et al. 1996, Natarajan et al. 1998, 
Geiger & Schneider 1998). For a given mass distribution, 
Kneib et al (1996) demonstrated that a redshift can be 
estimated if one can measure accurately the shape of an 
individual arclet. In order to validate this method and 
study its biases, extensive programs of gravitational arclet 
spectroscopy have been undertaken. In particular, Ebbels 
et al. (1998) have measured 19 redshifts of lensed objects 
behind Abell 2218. Most of them confirm the lens redshift 
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prediction, and allow the accuracy of the mass model to 
be improved. Similarly, Bezecourt and Soucail (1997) have 
started the spectroscopy of arclets in Abell 2390, which 
has been used to constrain the mass distribution in this 
cluster with a great accuracy (Kneib et al. 1998). 

Using these accurate cluster mass models and a spec- 
trophotometric description of galaxy evolution (Bruzual 
& Chariot 1993, Pozzetti et al. 1996), Bezecourt et 
al. (1998a) have predicted the expected arclet num- 
ber counts and their redshift distribution. This model 
presents many improvements with respect to previous 
work (e.g. Nemiroff & Dekel 1989, Wu & Hammer 1993, 
Grossman & Saha 1994, Hattori et al. 1997) as it includes 
many observational limits such as magnitude ranges, sur- 
face brightness cut-off or a choice of the optical waveband, 
and this for any mass distribution, regardless of its com- 
plexity. 

Abell 2218 is the first cluster where the number 
counts and redshift distribution of the background arclets 
have been examined in detail (Bezecourt et al. 1998a, 
Ebbels et al. 1998). We propose in this paper to further 
extend this study to another well known cluster lens, 
namely Abell 370. Its mass distribution was first accu- 
rately derived by Kneib et al. (1993) [hereafter K93] who 
showed that the mass model has to be bimodal in shape in 
order to accommodate the gravitational pair B2-B3. This 
was later confirmed by the HST/WFPC1 observation de- 
scribed in Smail et al. (1996) and the X-ray map of the 
cluster, displaying a bimodal shape compatible with the 
lens model (Fort & Mellier 1994). 

We present new observations of the cluster Abell 370 in 
Section 2: a deep HST/WFPC2 image and spectroscopic 
data. Section 3 discusses the new lensing mass model. In 
section 4, we use an improved version of the code devel- 
oped by Bezecourt et al. (1998a) to determine the ex- 
pected counts and redshift distribution of arclets in Abell 
370. Our analysis explores two different scenarios of galaxy 
evolution to study their differences, and compute the de- 
pletion curves of background number counts at different 
wavebands. Discussion, conclusions and further prospects 
to constrain galaxy evolution through lenses are presented 
in Section 5. Throughout the paper, we use a Hubble con- 
stant of H = 50kms _1 Mpc _1 , with £1 = or 1 and 
A = 0. 

2. New observational data 

2.1. HST observations and photometry 

Abell 370 was observed with the HST/WFPC2 camera 
with the F675W filter [ID: 6003, P.I.:R.P. Saglia], re- 
sulting in a reasonably deep image with a total expo- 
sure time of T exp = 5.6 ksec. These data were provided 
by the ESO/ST-ECF Science Archive Facility, Garching, 
Germany and were reduced with standard reduction pro- 
cedures. We used the IRAF/STSDAS packages for the im- 



age combinations after centering and cosmic ray cleaning. 
The absolute photometry was obtained using magnitude 
zero-points given in Holtzman et al. (1995). The photom- 
etry of the field was obtained with the Sextractor package 
(Bcrtin & Arnouts 1996). A criterion of 12 contiguous pix- 
els above the given detection threshold was used to iden- 
tify an object. The 1 a detection limit is i?675iv = 24.9 
mag arcscc~ 2 . From the magnitude histogram obtained 
from the catalog, we estimate the completeness limit to 
be R&75W — 25.5. 

We also built a sample of arclets for the purpose of 
our study of their photometric and statistical properties. 
To define an arclet we imposed the following selection cri- 
teria: at least 12 contiguous pixels above 2a of the local 
sky level, an axis ratio greater than 2, a central surface 
brightness lower than ^675W = 25.5 mag. arcsec -2 and 
a magnitude range 20 < Rqi^w < 26. The final sample 
contains 81 arclets and their magnitude histogram is dis- 
cussed in Section 4.2. 

2.2. Identification of multiple images and arclets 

Abell 370 (z=0.37) is a rich optical cluster dominated by 
two giant elliptical galaxies identified as #20 and #35, 
following the numbering of Mellier et al. (1988). A set 
of multiple image candidates and gravitational arcs are 
identified on the WFPC2/F675W image and are displayed 
in Figure 1 . Their photometric and geometrical properties 
are also summarized in Table 1. We now discuss them in 
detail: 

AO : Near galaxy #35 is located the spectacular gi- 
ant arc (z=0.724) initially detected by Soucail et al. 
(1987). From the high resolution WFPC2 image more 
details on its morphology are clearly visible, suggesting 
that it is a gravitationally lensed image of a spiral galaxy 
(Soucail et al. 1998a). 

Al to A6 : A first set of arclets (labelled Al to A6) was 
detected from ground based images by Fort et al. (1988) 
and further discussed by Kneib et al. (1994). Most of them 
are blue objects, but none have a spectroscopic redshift 
yet. A5 is the most extended one and presents very blue 
colors and a strong dimming in the reddest bands, sug- 
gesting a young star forming galaxy. Despite deep spectro- 
scopic data for A5 (Mellier et al. 1991) no significant emis- 
sion line has been identified suggesting that 1 < ZA5 < 2.2. 
Arclets A3 to A6 are probably single images of background 
lensed galaxies in view of their location in the lens plane. 
Arclets Al and A2 may be multiple images of the same 
source. 

The multiple images B, C and D : K93 demonstrated 
that the B2-B3 objects (Fort et al. 1988) correspond to 
a gravitational pair with a counter image labeled as B4. 
The B2-B3 pair and the arc AO were used to constrain 
the K93 model which showed a bimodal mass distribution. 
This model proposed that objects C1-C2-C3 and D1-D2, 
identified on high quality ground-based images, were also 
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Fig. 1. Detailed view of the multiple image candidates detected in the WFPC2/F675W image. B2-B3-B4 is a triple image 
configuration, as well as Cf-C2-C4. R1/R2 is a radial arc. E1-E2-E3 is also a triple configuration with a clear inversion of 
parity between E2 and E1/E3 (see text for more details). 



multiple image systems. A redshift estimate based on the 
mass model was proposed for each of the pairs: zb — 0.865, 
zc = 0.81 and zd = 0.95. The reality of B4 was confirmed 
a posteriori by the HST/WFPC1 data (Smail ct al. 1996). 
C3 is likely to be a wrong identification of the counter 
image of C1-C2, and we denote C4 the correct counter 
image used in the present model. 

The radial arc R : In the HST/WFPC1 image, Smail 
et al. (1996) discovered a radial arc candidate R (R1-R2) 
located close to galaxy #35. They modeled it as a 5-image 
configuration, and predicted a redshift of z ~ 1.3 using the 
K93 model. This radial arc is well identified in the F675W 
image and is clearly the merging of 2 images. 
New multiple images : The detailed insight of the mass 
model and the exquisite HST resolution allow the identifi- 
cation of other multiple image candidates. We only discuss 
here the E1-E2-E3 multiple configuration as it presents a 
characteristic inversion of parity as expected from lensing 
theory (see also Smail et al. 1995 and Colley et al. 1996 
for other spectacular examples). Other multiple image 
candidates will be discussed in a forthcoming paper. 

2.3. Spectroscopic Observations 

Spectroscopic data were acquired at the 3.6m telescope of 
La Silla (ESO) with EFOSC on October 1996. A long slit 
of 1.5" width was positioned along the two objects B2 and 
B3 for a total integration time of 3 hours with an average 
seeing of 1.9". Following the predictions of K93, the [On] 
emission line was expected to lie at A ~ 6950A. Thus the 
R300 grism was used, providing a useful spectral range 
6000A-9000A and a dispersion of 7.5A/pixel. The data 
were reduced with standard procedures for flat fielding, 
wavelength and flux calibration in the IRAF environment. 



Table 1. Main geometrical and photometric properties of the 
multiple images candidates identified on the F675W HST im- 
age. The origin of the coordinates is the center of galaxy #35 
and the XY orientation is the CCD one, i.e. North is to the 
top with an angle of 85° 41' clockwise from the X axis and 
East is to the left. Coordinates are in arcseconds. e is the el- 
lipticity of the objects {i.e. 1 - b/a where b/a is the axis ratio 
of the isophotes) and 9 is the orientation of the isophote with 
respect to the X axis. Surface photometry was computed on 
the two brightest elliptical galaxy using the ellipse package in 
STSDAS. The two galaxies are well fitted by a de Vaucouleurs 
law with R e = Z2.&h^ kpc for #20 and R e = Uh^ kpc for 
#35. 



Object 


X(") 


Y(") 


e 


6 


R675W 


Al 


16.35 


57.30 


0.74 


-23.8 


24.83 


A2 


5.56 


60.43 


0.75 


-4.6 


24.38 


B2 


8.46 


21.14 


0.76 


-0.2 


23.08 


B3 


13.80 


20.58 


0.74 


-15.7 


23.23 


B4 


-19.45 


20.98 


0.42 


-41.6 


23.82 


CI 


6.46 


6.50 


0.17 


4.6 


24.05 


C2 


10.19 


4.96 


0.10 


63.2 


24.00 


C3 


-21.82 


2.78 


0.14 


-60.9 


23.99 


R1/R2 


5.27 


7.03 


0.76 


76.4 


23.98 


El 


32.86 


18.62 


0.51 


74.6 


24.83 


E2 


-31.22 


19.15 


0.68 


64.4 


24.53 


E3 


0.54 


22.21 


0.66 


-69.8 


24.44 


#20 
# 35 


3.09 
0.00 


37.85 
0.00 


0.36 
0.45 


-70±3 
-77±4 


16.9 
17.1 
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Sky subtraction was performed on the 2D image, and the 2.4- A closer look at the B2-B3-B4 multiple image 

final spectra were extracted with an optimal extraction 

algorithm. 




3500 4000 4500 5000 




6000 6500 7000 7500 8000 8500 9000 
wavelength (A) 



Fig. 3. Spectra of objects B2 and B3 co-added and flux cali- 
brated. The ordinate is F\ in arbitrary units with zero level at 
the bottom of the graph. The wavelength scale is the observed 
one at the bottom, and the rest frame one at the top. The two 
emission lines detected in the spectra are quoted with their 
identification, as well as the location of the absorption H and 
K lines, not clearly detected. 

Both objects B2 and B3 show a prominent emission 
line at A = 6728A and A = 6727A respectively on the two 
dimensional spectra, which we identify with [On] 3727A 
(Figure 2 and 3). Another feature appears at A = 9045A 
and A = 9042A respectively, identified with [Om] 5007A 
thus confirming the redshifts. We find zb2 = 0.8058 and 
zb3 = 0.8054, clearly demonstrating the similarity of the 
two objects within the error bars, and in good agreement 
with the z = 0.865 prediction. This definitely confirms the 
nature of these objects as a multiply imaged galaxy. 

Although, arclets Al and A2 were observed during 4.5 
hours, no continuum nor emission line could be detected. 



Fig. 4. Multicolor photometry of object B3 in U,B,R,I,J 
and K' . The best fit of the data points corresponds to a star 
forming galaxy with a constant star formation rate seen at 
t — 1.7 Gyr, with Z = Zq and a null absorption. Magni- 
tudes are: U = 22.43 ± 0.1, B = 24.58 ± 0.1, R = 23.84 ± 0.2, 
J = 22.08 ± 0.2, K' = 21.26 ± 0.5. 

This multiple image is one of the bluest arclets ob- 
served in this field. Accurate photometry is available in 
many filters extending from U to K' with data coming 
from a U HST image with the F336W filter (ID: 5709, 
P.I. J.M Deharveng, Bezecourt et al. 1998b, hereafter Pa- 
per II), B, R and I images from CFHT (Kncib et al. 1994) 
and unpublished CFHT infrared images in J and K' bands 
(Le Borgne and Picat, private communication). B2 lies too 
close to galaxy #20 and is contaminated by its envelope 
at long wavelengths, hence we only concentrate on B3. 

We compute the spectral energy distribution (SED) 
of B3, corrected from its redshift, in order to study the 
stellar content of the galaxy. We compare and fit the B3 
SED with synthetic ones from the Bruzual and Chariot 
spectrophotometric evolutionary code (GISSEL96). This 
code can provide SEDs for different metallicities ranging 
from Z Q /50 to 5 Zq, different spectroscopic galaxy type 
determined by the star formation history (burst, ellipti- 
cal, spiral and irregular), and different internal absorp- 
tion by dust can be studied. The extinction law used is 
the one proposed by Calzetti (1997) which is similar to 
what is observed in our galaxy without the 2175A bump. 
We explore the 3-dimensional parameter space (spectro- 
scopic type, metallicity and absorption, E(B — V) from 
to 0.3 mag.) to fit the B3 photometry by the model SEDs. 

The main constraint comes from the high UV flux 
which requires a weak internal absorption (E(B — V) < 
0.2mag) and a low metallicity (Z < Zq) (Figure 4). De- 
spite this, the fit does not reproduce the UV part very 
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well. A possibility is that we are seeing the superposition 
of a very recent burst of star formation on an older stellar 
population. A strong starburst alone would show a much 
fainter flux at red and infrared wavelengths, hence the ob- 
served IR emission must come from an already existing 
old population. The B multiple image is similar to arclet 
A5 which is also very bright in the UV but much fainter 
in the IR (R — K' = 2.6 for B3 and R - K' < 1.9 for A5). 
More precisely, assuming a redshift of z—1.3 (Mellier ct 
al. 1991), we have similarly fitted the A5 SED. It shows 
evidences for a low absorption (E(B — V) = 0.1 mag.) and 
a very low mctallicity (Z = Z Q /50) (Figure 5). The fact 
that the best-fit metallicity (Zq /50) is the lower bound al- 
lowed by the code should not be a problem as a satisfying 
fit is also obtained for Z & /5 with comparable significance. 
The trend towards low Z is then stable and not too fast. 
The corresponding age is 0.3 Gyr, drawing the portrait of 
a very young object in an active phase of star formation. 
An additional old population is not required for A5 which 
is the major difference with object B. 
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Fig. 5. Multicolor photometry of object A5 in U, B, R, I, and J 
(undetected in K'). The best fit of the data points corresponds 
to a burst of star formation redshifted at z = 1.3 and seen at 
t = 0.3 Gyr, with Z = Z Q /50 and E(B-V) = 0.1 mag. Magni- 
tudes are: U = 21.70±0.15, B = 22.90±0.04, R = 22.26±0.06, 
I = 21.67 ± 0.09, J = 20.07 ± 0.11. 



These two objects B3 and A5 are typical examples of 
the population revealed by UV imaging. The UV selec- 
tion exhibits young objects with low mctallicity and low 
absorption in the redshift range [0.5,2.0]. 

Such detailed stellar population study, could in prin- 
ciple be applied on any other arclets in the field, for 
which redshifts have been spectroscopically confirmed or 
are strongly constrained by lensing and photometric red- 
shift techniques (Pello et al. 1998a). 
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3. An improved lens model 

Using both the giant arc AO and the B2-B3-B4 triple sys- 
tem, K93 showed that a bimodal mass distribution was an 
adequate fit to these multiple image constraints (see also 
AbdelSalam et al. 1998 for a different approach). Here, 
we improve this gross picture by taking into account the 
contribution of the cluster galaxies (mainly E/S0 galax- 
ies in the cluster core) in a similar way as Kneib et al. 
(1996) and Natarajan et al. (1998). All mass components 
(galaxies+clusters) have been assumed to follow a trun- 
cated pseudo isothermal mass distribution (PIEMD, Kas- 
siola and Kovncr 1993, Hjorth & Kneib 1998). 

Following Kneib et al. (1996), for each galaxy halo, the 
velocity dispersion a , the truncation radius r t and the 
core radius tq are scaled to the galaxy luminosity com- 
puted from the observed F675W magnitude. The scaling 
relations used for the galaxy halos are: 

re = r,(-j , (2) 

ro = r »(£) 2 . (3) 

The scaling relations adopted are motivated by the prop- 
erties of the fundamental plane (FP) and are similar to 
the one used by Brainerd et al. (1996). In particular, the 
exponent 0.8 in eq. (2) leads to a total mass-to-light ra- 
tio that scales with L 3 in agreement with the observed 
correlation of the FP (e.g. Jorgensen, Franx & Kjaergaard 
1996). Brighter galaxies have more extended and therefore 
more massive dark halos. 

The orientation and ellipticity of the galaxy halos 
are taken from the observed values of the light distri- 
bution while (To* an d ft* are optimized, ro* is fixed at 
0.15 kpc. To model the cluster components we considered 
two 'large scale' mass distributions (represented by 2 trun- 
cated PIEMD mass distribution) centered on galaxies #20 
and #35. Their orientations, ellipticities, velocity disper- 
sions and core radius are left as free parameters. 

The constraints imposed in the optimization procedure 
are the multiple images described in section 2.2, along 
with the redshifts of AO and of the B2-B3-B4 triplet. The 
best fit parameters are summarized in Table 2 (x 2 =4.5). 
We also derive a "lensing" redshift estimate for the mul- 
tiple images. The C system is at zq — 0.75 ± 0.1, D is at 
zd = 0.85 ± 0.1 and E is at ze — 1.3 ± 0.1. R is moved to 
larger redshift with zr = 1.7 ±0.2 and A1/A2 could form 
a gravitational pair at zai/A2 = 1-4 ± 0.2. Wc computed 
the total projected mass in different apertures centered 
at the barycentre of the #20 and #35 galaxies. We find 
that within 75, 150 and 300 kpc the total mass is, respec- 
tively, for this modeling and the previous one (K93): 0.5 




Fig. 6. Central part of Abell 370 seen with WFPC2 in F675W. Overlaid is the mass distribution (black thin lines) and the 
critical lines at z — 0.806 (thick white lines) and z — 4 (thick black lines). The most important arclets and multiple images are 
shown. 



(0.45) ±0.05 10 14 M , 1.8 (1.6) ±0.1 10 14 M and 4.8 
(4.3) ±0.15 10 14 M . Furthermore, 5% of the total mass 
is retained in cluster galaxy halos. The total mass-to- light 
ratio is -180 (M/L v ) (out to 300kpc, 16O(M/L y ) for 
the K93 modeling), and —9 (M/Ly) for a L* galaxy ha- 
los. These results are similar to the one found in other 
cluster lenses like Abell 2218, 2390 and AC114 (Kneib et 
al 1996, 1998; Natarajan et al 1998). 



4. Statistical properties of faint galaxies 

4.1. Modeling the number counts of tensed galaxies 

The number counts of gravitational arclets in the field of a 
massive cluster of galaxies are a competition between the 
magnification of the luminosity by the cluster potential 
that makes more objects visible and the surface dilution 
that decreases the surface density of arclets by the same 
factor as for the magnification. 
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Table 2. Model parameters for A370 potential 
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the complex relation given by Kneib et al. (1996) that re- 
lates the image deformation tj as a function of the source 
deformation Ts and the strength of the potential T pot : 



sgn(det A 1 )t i = t s + T pot (S s + SR(t s9* po t)) 



The number of arcs brighter than magnitude m with 
an axis ratio greater than q m in and a surface brigthness 
brighter than /iq within the field of a cluster of galaxies is: 



N(m,q min ,fj, ) 



oo 



S(q,z) Jl 



$ i (L,z)dLdq^dz 



(4) 



The sum is over the different morphological types i. zi 
is the lens redshift and z max ((io,i) is the redshift cutoff 
corresponding to the limit in central surface brightness /iq. 
S(q, z) is the angular area in the source plane (at redshift 
z) that gives arcs with an axis ratio between q and q + dq. 
$>i(L,z) is the luminosity function at redshift z for each 
morphological type. 

It was required as a preliminary step that counts in 
wavebands from U to K and redshift distributions in B 
and / in empty field are correctly reproduced. We used the 
model for galaxy evolution of Bruzual and Chariot (1993) 
with the prescriptions of Pozzetti et al. (1996) for the op- 
timisation of the parameters representing the field galaxy 
distributions (see Bezecourt et al. 1998a for details). 

The weak point in the Bezecourt et al. 1998a model, 
is the approximation of circular sources used for the cal- 
culation of arclets axis ratio. Indeed, this may underesti- 
mate the number of galaxies with axis ratio larger than 
a certain threshold especially at low redshifts and then 
bias the predicted redshift distribution. To provide a more 
accurate and reliable redshift distribution, we now take 
into account the galaxy ellipticity distribution as given by 
Ebbcls (1998) (this distribution is based on the analysis 
of HST/MDS-like images, and is therefore appropriate to 
our study): 



p(r) oc t exp 



where r = 



(o.o3e) 



0.54 



(5) 



2ab 



for an elliptical object with semi ma- 
jor axis a and semi minor axis b. The complex deformation 
t is defined by: r = r cos(2#) + it sin(20) where 9 is the 
orientation of the major axis. At each point in the image 
plane (and at each redshift step for the sources), we es- 
timate the fraction of the background galaxies which can 
give a lensed image more distorted than a fixed limit. We 
first consider the effect of the lens deformation through 



(6) 



where 8 = v(T+ t 2 ) = 1 + gr* (see Kneib et al 1996 for 
a full description of this formalism). The term 5ft(rs<7* ot ) 
is a correction for strong lensing only. A lower threshold 
in the observed axis ratio (a/6 > 2 here) corresponds to 
a lower threshold in the deformation (rj > 0.75) what- 
ever the position angle of the image. Within this limit, we 
then scan all the allowed solutions for ts and the source 
position angle 9s- The fraction F(xi,yi, z) of galaxies at 
the location (xi,yi) in the image plane and at redshift z 
which fill the condition is computed using the probability 
distribution of equation (4) for rj and an average over 
the position angle 9s- 



T{xi,yi,z) = / / p(T S )dT S d9 s 

Z7T JO Jts I r/>0.75 



(7) 



4-2. Comparison between the observed and predicted ar- 
clets number counts 




F675W 

Fig. 7. Number counts of arclets in A370 in the F675W 
HST image with the following selection criteria: a/b > 2 and 
M675W < 25.5. Observed counts (o) are compared to two num- 
ber counts models with respectively qo = 0.5 (solid line) and 
go = (dashed line). Error bars are statistical poissonian un- 
certainties. 



The magnitude histogram of the arclets observed 
in A370 is shown in Figure 7. As already stressed in 
Bezecourt et al. (1998a), we find that the addition of 
galaxy scale components and the effect of the source ellip- 
ticity distribution significantly increase arclets counts by 
a factor of 1.4 (qo = 0.5) and 1.7 (qo = 0) with respect 
to a K93 model with the assumption of circular sources. 
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Such increase in the arclets counts was already observed in 
Bezecourt ct al. (1998a) when using a detailed mass distri- 
bution but only with circular sources. Considering an cllip- 
ticity distribution for the galaxies also increases the arclet 
counts and strongly modifies the redshift distribution as 
discussed below. With this refinement, the arclet counts 
prediction is now more consistent with the data over the 
whole magnitude range. Clearly detailed lens models and 
proper assumptions on the source ellipticity distribution 
are mandatory to explore arclets counts. 



4-3. Redshift distribution of arclets 



In order to sec how stable the redshift distribution is 
with respect to the ellipticity distribution, another p(r) is 
considered following Ebbels (1998): 



p(t) oc t exp 




Fig. 8. Redshift distribution of arclets in A370 for the follow- 
ing selection criteria: R&nw < 23.5, /j,r < 25.5 and a/b > 2. 
The solid line corresponds to the model qo = 0.5 and the 
dashed line is for qo = 0. The dotted line corresponds to the 
ellipticity distribution of equation 8 with q = 0.5. 



The next step is to consider the redshift distribution 
of the arclets. Figure 8 shows the predicted one with 
the selection criteria close to the observational limits for 
faint object spectroscopy (i?675 < 23.5, < 25.5 and 
a/b > 2). The comparison with real data is not easy at 
present because no well defined sample of arclets, complete 
in magnitude, has yet been spectroscopically explored. 
The redshift distributions displayed in Figure 8 present 
a prominent peak at z ~ 0.6 and a secondary one at z > 2 
due essentially to the contribution of bright and young 
Elliptical-type galaxies in their process of strong star for- 
mation. This second peak is strongly attenuated com- 
pared to what was discussed in Bezecourt et al. (1998a) 
where there was clearly a strong excess in the high redshift 
tail. Taking into account the ellipticity distribution of the 
sources has fixed this problem. 



(o.2o) 



0.8F, 



(8) 



representative of objects with 18 < / < 25.5 while equa- 
tion 5 was derived with bright objects (18 < I < 22). As 
this distribution corresponds to more elongated objects 
than in equation 5, more arclets appear at low z (Figure 
8, dotted line). However, the total number of arclets is 
increased by only 11%. 

Analysis of Figure 8 suggests several comments: 

— Compared with the redshift distribution of field galax- 
ies within the same observational limits, the arclet dis- 
tribution is biased towards more distant objects. We 
can take advantage of this modification of the redshift 
distribution to select more distant galaxies in a redshift 
survey sample. This is particularly true near the cen- 
tral part of the lens where very high redshift galaxies 
(z > 2.5) have already been found (Ebbels et al. 1996, 
Trager et al. 1997, Franx et al. 1997, 
Frye & Broadhurst 1998, Pello et al. 1998a). 

— The two galaxy evolution models presented in this pa- 
per (qo = and qo = 0.5) have been validated in empty 
fields where the number counts and redshift distribu- 
tion of galaxies were compatible with the observed 
ones. However they present significant differences at 
high redshift for arclets. The q = model with "Pure 
Luminosity Evolution" (PLE) predicts more high red- 
shift arclets than the qo = 0.5 model where number 
density evolution is included. This different behavior 
of the redshift distribution at z > 2 is encouraging 
as it could be a way to distinguish the two scenarios 
by analyzing the redshift distribution of arclets with 
z > 2. Although it is presently difficult to speculate on 
the true fraction of very high redshift objects in well 
defined samples of arclets, this population of very high 
redshift arclets does not seem to dominate the actual 
samples. Hence we favor the qo = 0.5 model corre- 
sponding to a merging scenario of galaxy formation 
and evolution. 

— Because the high redshift domain is sensitive to many 
uncertainties included in the evolutionary code (dust 
obscuration in the UV, uncertainties in the slope of 
the IMF for massive stars or in their UV tracks, influ- 
ence of the short time scale phenomena), we must be 
cautious on the above conclusions. Furthermore, the 
galaxy evolution models assume that all the galaxies 
form instantaneously at a given redshift, which is a 
clear limitation of these models (we have Zf olm = 4.5 
for qo = 0, which gives for the present day galaxies an 
age of 16 Gyr, and Z[ orm = 5 for qo = 0.5, correspond- 
ing to an age of 12.2 Gyr). A more phenomenological 
model of galaxy formation that follows in detail the 
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mass evolution of galaxy sub units (Baugh et al. 1998) 
may improve this simple description, and would likely 
be the best to be compared with the current data. 

4-4- Depletion curves in Abell 370 




r (") 

Fig. 9. Depletion of lensed object in different filters and lim- 
iting magnitudes: U < 23 (-•••-), B < 24.5 ( — ), R < 24 

( ), / < 23 (- • -) and K < 21 (• • •) for go = 0. The 

small excess near r = 40" is due to the local magnification 
enhancement of an individual cluster galaxy. 

For a simple description of the background popula- 
tion, it is easy to show that the number density of objects 
brighter than a given magnitude m behind a lens is: 

N(< m,A) = N (< m)A 2 - 5a - 1 (9) 

where No(< m) is the number density in blank field, A is 
the magnification and a is the logarithmic slope of field 
number counts. An excess or a lack of objects is then ex- 
pected for a slope steeper or shallower than 0.4 respec- 
tively. In practice, optical number counts of field galax- 
ies show a slope smaller than 0.4 in nearly all wavebands 
(except in U and B at relatively bright magnitudes), so 
a "depletion" is expected in most cases, being more pro- 
nounced at longer wavelengths (Fort et al. 1996, Taylor et 
al. 1998). 

With our model, it is possible to compute radial de- 
pletion curves instead of global number counts, for any 
filter or magnitude range. A clear illustration of the wave- 
length dependence of the magnification bias is given by 
the ratio of the number of lensed objects expected in a 
given area in the field of Abell 370 over the number of 
field galaxies in the same area (Figure 9). One can note 
that the predicted intensity of the depletion is higher at 
longer wavelengths, as expected, because of the shallower 
slope of the field galaxies counts. Moreover, because of the 
flattening of the counts at faint magnitudes, the depletion 



is also very sensitive to the magnitude threshold (see the I 
curve in Figure 9). UV counts have a slope larger than 0.4 
at bright magnitudes (U < 23) but it quickly decreases 
at fainter levels. This fast flattening is due to the Lyman 
break that goes through the red limit of the U filter at 
redshift z ~ 2.5. Hence, the lack of objects at faint magni- 
tudes in U produces the depletion curve shown in Figure 
9 (see Paper II for more details of the UV modeling). 

The detection of this magnification bias is highly de- 
pendant on the poissonian noise of the background sources 
and the contamination by cluster members and foreground 
objects. This is critical in cluster cores where the den- 
sity of objects is very low, of the order of a few units per 
arcmin 2 . Hence, these very poor statistics cannot bring 
valuable information on the sources redshift distribution 
as most of the information is expected in the dip of the de- 
pletion curve. Only massive clusters are able to show such 
a depletion curve but the recovery of the sources redshift 
distribution seems out of reach with present datasets. 

5. Conclusions 

In this paper, we have presented a new analysis of the 
cluster lens Abell 370. Thanks to the measurement of 
the redshift for the B2-B3 gravitational pair (z = 0.806) 
and the identification of several new multiple arclets in 
the WFPC2 image, a more accurate and well constrained 
model of the mass distribution in Abell 370 is proposed, 
including galaxy scale components. We studied the spec- 
tral energy distribution (SED) of arclets B3 and A5 which 
are found to be both young, low metallicity star forming 
objects without strong interstellar extinction. 

The lens model has been used to study the background 
population of galaxies. Taking the galaxy ellipticity dis- 
tribution into account induces significant changes in the 
predicted redshift distribution of arclets. The excess of 
very distant sources found in previous analyses is now 
strongly attenuated and arclets at low redshift are recov- 
ered. For the two prescriptions used on galaxy evolution, 
although they both reproduce well the number counts in 
empty fields and in cluster lenses, the "Pure Luminosity 
Evolution" model in a low density universe over predicts 
the number of high redshift (z > 2) galaxies, compared to 
what is currently observed. If this effect is observationally 
confirmed it may constrain the number density evolution 
of galaxies in order to interpret the redshift distribution 
of arclets. 

To understand further the properties of the back- 
ground population through cluster lenses, deep dedicated 
cluster surveys are needed in order to enlarge the num- 
ber of arclets and have significant statistical properties. A 
multi-color approach is favored as it can help in constrain- 
ing the redshift distribution. This is particularly true for 
those unresolved sources for which no shape information is 
available and only the spectral energy distribution can dis- 
criminate between background objects and cluster mem- 
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bers. In addition, detailed study of their spectral energy 
distribution will be useful to characterize properties like 
the luminosity function, dust extinction and star forma- 
tion histories of distant galaxies (e.g. Pello et al. 1998b). 
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